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В данной работе рассматривается задача применения метода вычислительного эксперимента к 
исследованию процесса сушки сыпучих материалов в «кипящем слое» как вероятностного явления. 
На основе использования математических моделей диффузионно-переносного движения, нагрева и 
обезвоживания частиц исследованы вероятностные характеристики времени пребывания, температур 
и степени обезвоживания материала в «псевдоожиженном» слое. Полученные результаты показывают 
преимущества способа по сравнению с другими технологиями. 
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Введение 


Особенностью технологического процесса обработки материалов в псевдоожижен- 
ном слое является случайная природа времени пребывания частиц в рабочей зоне, 
обусловленная хаотическими движениями и соударением частиц. 

В работах В.Н. Ткаченко [1-4], обобщающих теоретические результаты ряда авторов, 
в частности В.В. Кафарова [5], Г. Крамера [6] и других, предложен метод вычислитель- 
ного эксперимента для исследования стохастических параметров процесса воздействия на 
материалы в кипящем слое. Применение данного метода дает дополнительные возмож- 
ности идентификации математических моделей процессов обезвоживания сыпучих мате- 
риалов в псевдоожиженном слое. Вычисление показателей плотности распределения 
частиц по времени пребывания в псевдоожиженном слое. 

Для моделирования процесса изменения концентрации частиц в псевдоожи- 
женном слое используется математическая модель диффузии частиц в конвективном 
потоке [1], [5]. Уравнения моделируемого процесса имеют вид: 


2 
ас _ 58'С_ 8С 


ОЕ дх” (6.4 (1) 
р 9С(0,т) ая 
ОХ (2) 
ОС(Г,т) о 
ОХ (3) 
С(0,Х) = С, (Х), (4) 


где 

С(х,х) — безразмерная концентрация частиц в точке х, в момент времени т; 

Р — коэффициент диффузии частиц в слое; 

у — расходная скорость потока; 

С0(Х) - начальная концентрация частиц. 

Для реализации математической модели (1) — (4) на ЭВМ могут быть использо- 
ваны конечно-разностные методы с явной или неявной схемой аппроксимации и их 
модификации [1], [7-9]. С помощью метода расщепления модель (1) — (4) может быть 
представлена элементарными физическими процессами (диффузией и переносом 
массы), для расчета которых накоплен значительный опыт. Простейшей для реализации 
на ЭВМ является явная схема, в соответствии с которой краевая задача представ- 
ляется в виде следующей конечно-разностной: 


й т т т 7-1 й 7-1 
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Шаг по времени Вт выбирается из условия устойчивости, которое получено по 
основе анализа уравнения (5) и имеет следующий вид: 
2 
Й, < о 
2О+уйу 

Конечно-разностное уравнение (5) аппроксимирует исходное уравнение по схеме 
односторонней разности против потока с первым порядком точности по переменным Х.л. 
Для определения параметров конечно-разностной схемы (Вх ‚ Вт) отладки алгоритмов 
расчета краевой задачи удобно пользоваться интегральными характеристиками про- 
цесса, основанными на законах сохранения массы. 

При равномерно распределенной загрузке материала для больших коэффи- 
циентов диффузии (О > 0,01) режим движения частиц соответствует идеальному 
перемешиванию, а закон распределения частиц по времени пребывания близок к 
экспоненциальному. При малых же коэффициентах диффузии (О < 0,01) плотность 
распределения частиц превращается в функцию, напоминающую агс®, причем, чем 
меньше О, тем длиннее площадка равномерной плотности и распределение стремится к 
прямоугольному. Проведены исследования вероятностных характеристик времени 
пребывания в диапазоне значений коэффициентов диффузии, полученных опытным 
путем для различных материалов, обрабатываемых в псевдоожиженном слое [1], [5]. 
В табл. | приведены результаты численного исследования вероятностных характе- 
ристик времени пребывания частиц в рабочей камере при различных значениях 
параметров процесса. 


Таблица 1 — Зависимость вероятностных характеристик времени пребывания 
частиц от величины коэффициента диффузии 


1 
3 


0 
1 
9 
0 
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В ячейках таблицы первое значение обозначает математическое ожидание, второе — 
среднеквадратичное отклонение, третье — коэффициент вариации. 

В результате исследования установлено, что математическое ожидание времени 
пребывания частицы в активной зоне составляет более 11 мин, в переходной зоне — 
более 13 мин, что подтверждает как результаты моделирования, так и положение о 
том, что сушка в «кипящем слое» обеспечивает более длительный контакт материала 
с сушильным агентом по сравнению с барабанными сушилками, трубами-сушилка- 
ми и другими схемами, и по этой причине является наиболее эффективной современ- 
ной технологией обезвоживания. 
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Исследование температуры и степени обезвоживания частиц в псевдоожижен- 
ном слое как стохастических величин. 

В простейшей одномерной постановке основные процессы теплообмена описывает- 
ся следующей системой обыкновенных дифференциальных уравнений: 


а-эс тб Тат т) 
У, тел м. 2 Су 5 5 


ат, 
ОХ ОХ (9) 
Эт КУР 
вы 3 ы Е 
ВР Г.) р (Г, Ти) 
0<Х <[,0<1<&, 


(10) 
ОТ, (7,0) _ 0 
р на еа. 

5 ох 5 25 5 5 (11) 

ета, 
Эа 5 575 5 (12) 

2150 

Т(Р)=То. ие 


Характерной особенностью математического моделирования теплообмена в 
противотоке является тот факт, что известными оказываются, как правило, тем- 
пературы материала и газа на входе, задаваемые на разных концах зоны теплообмена. 
При численном интегрировании системы (9) — (13), начиная в точке х = 0, либо в 
обратном направлении при х = Г, оказываются неизвестными Т5(0), либо Т$. Целью 
математического моделирования работы реактора является расчет функции Т5(Х), 
То(Х), для Х: 0 <Х <Г и значения температуры материала на выходе — для заданных 
расходных скоростей у$, уз. Для расчета требуемых величин необходимо решить 
некоторую краевую задачу, для чего дополним систему уравнений (9) — (13) началь- 
ным условием для температуры газа Т(0) в точке х = 0, таким, чтобы при интегри- 
ровании системы получить равенство рассчитанной в точке х = Г. температуры газа 
Тег (Г.) заданной начальным условием (13): 

г 0 
ИН 

Определение начального условия Т°(0), для которого выполняется равенство 
(13), сводится к нахождению корня некоторого уравнения, заданного алгоритмически. С 
этой целью система уравнений аппроксимируется какой-либо конечно-разностной 
схемой, обеспечивающей достаточно высокую точность расчета искомых функций 
[8], после чего решение краевой задачи может быть получено одним из численных 
методов поиска корней алгебраических и трансцендентных уравнений (методы Ньютона, 
хорд, деление отрезка пополам и др.). 

Расчет вероятностных характеристик осуществляется в едином алгоритме интегри- 
рования краевой задачи, на каждом шаге которого накапливаются суммы требуемых 
величин, умноженных на плотность распределения. 

На основе законов распределения частиц по времени пребывания, полученных 
с помощью диффузионной модели, исследованы вероятностные характеристики 
температур и степени обезвоживания частиц в псевдоожиженном слое. Представляет 
практический интерес установить, в каких пределах среднего квадратического от- 
клонения времени пребывания, а значит, и при каких значениях коэффициента 
диффузии, равенство М[Т(т)]=Т(т,) все же можно считать приближенно выпол- 


НЯЮЩИМСЯ. 
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В табл. 2 приведены вероятностные характеристики температур и степени обезво- 
живания частиц в диапазоне значений коэффициентов диффузии, полученных опытным 
путем для различных материалов, обрабатываемых в псевдоожиженном слое [5]. 
С увеличением коэффициента диффузии растет среднее квадратическое отклонение 
(с.к.о.) времени пребывания при неизменном его математическом ожидании (тт = 40 мин). 
При этом изменяются не только дисперсии, но и математические ожидания всех характе- 
ристик частицы, что подтверждает нелинейность соответствующих функциональных 
зависимостей. Степень нелинейности различных характеристик частицы существен- 
но отличается. 


Таблица 2 — Вероятностные характеристики температур и степени 
обезвоживания частиц при различных значениях коэффициента диффузии и законах 
распределения времени пребывания 


С.к.о. Температура Среднемассовая Степень 
времени поверхности температура обезвоживания час- 
пребывания частицы ТИЦЫ 
(мин) 


Коэффициент 


диффузии 
м2 /с 


Аппроксимация нормальным законом распределения 


Режим идеального смещения (экспоненциальный закон) 


40,0 


Примечание. В таблице над разделительной чертой указана величина 
математического ожидания, под разделительной чертой — среднее квадратическое отклонение 
случайной величины. 


Максимальной нелинейностью обладает степень обезвоживания частицы. При 
О < 0,0001 (соответственно от < 7 минут) вероятностный расчет с допустимой по- 
грешностью можно заменить детерминированным расчетом температур и степени 
обезвоживания. При О > 0,0003 (от < 10 минут) игнорирование случайного характера 
времени пребывания частиц приводит уже к существенным погрешностям в числен- 
ной оценке степени обезвоживания частицы. 
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Проведен вычислительный эксперимент с аппроксимацией плотности распре- 
деления времени пребывания, полученной при О = 0,0003 (строка 4 табл. 2), нор- 
мальной плотностью для одинаковых значений математического ожидания и дисперсии. 

Сравнение результатов расчета вероятностных характеристик температур и 
степени обезвоживания частиц позволяет сделать вывод, что в классе унимодальных 
законов распределения частиц по времени пребывания математические ожидания 
характеристик частиц зависят от математического ожидания и дисперсии времени 
пребывания частиц и практически не зависят от моментов более высокого порядка, 
тогда как дисперсии зависят еще и от асимметрии и эксцесса. Практически это 
означает, то при оценке степени обезвоживания по заданным математическому ожида- 
нию и дисперсии времени пребывания частиц и неизвестной плотности распреде- 
ления, последнюю можно считать нормальной. 


Выводы 


Анализ вероятностных характеристик показывает, что при движении материала 
из зоны в зону дисперсия времени пребывания частиц нарастает, однако медленнее, 
чем его математическое ожидание, в связи с чем коэффициент вариации времени 
пребывания уменьшается, при этом закон распределения стремится к нормальному 
(асимметрия и эксцесс уменьшаются с увеличением числа пройденных зон). 

Сравнение результатов расчета вероятностных характеристик температур и степени 
обезвоживания частиц позволяет сделать вывод, что в классе унимодальных законов 
распределения частиц по времени пребывания математические ожидания характе- 
ристик частиц зависят от математического ожидания и дисперсии времени пребы- 
вания частиц и практически не зависят от моментов более высокого порядка, и это 
позволяет при оценке степени обезвоживания по заданным математическому ожиданию 
и дисперсии времени пребывания частиц плотность распределения можно считать 
нормальной. 
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Тре пуезиеайоп о} Капаот Сйагасчету ор ОшсЁзапа Маета[5 Огуия 


Ргосезб т «ВоШие Гауег» Бу Метоа о} Сотршег Ехрентет 

ВасКогоипа: ш отуеп \уогК Фе ргоМет оЁ аррИсаНоп о сотршег ехрегипепЕ Юг 
шуезИрайоп оЁР аисКзап@ таепа!5 Чгуше ргосез$ ш «БоШие 1ауег» аз гапдот 
рВепотепоп 15 сопз14егед. Оп Фе Базе оЁ аррИсаНоп оЁР тафетайса! то4е|$ оЁ т]ес{- 
фгапзРег тоуше, Беайпе ап ужег-ехфасипе Фе гапдот рагатеег$ ог ех1$Ипе Ите, 
{етрегабте ап Чгушо 4еотее оР таепа| ш рзечдо-Паш4 1ауег \аз шуезиоаеа. Те 
оМаштед гези 5 зВо\/ Фе ргеРегепсе оР тефо@$ сотрагей у апоег {есбпо]оз1е$. 

Маена5 ап ше#о4$: Те шуезИсайоп 15 деуо{е4 10 зоаЯоп оЁР а ргоет оЁ 
ре есйоп оЁ тео4$ оЁ еогейса! шуезИсайоп ап4 отоипше оё рагатаег$ оф ргосе$з 
ОЕ агуше оЁ роиг таепа]$ изше тефо4 оЁ "БоШие 1ауег", ргоу14е4 еНесйуепез$ оЁ 
фесбпо]оэу ап4 диаШу оф рго]есйпе. Ргеземе4 фе пе\у тафетайса| то4е|5 оЁ д&егиттеа 
буре, ехрегитеша! шуезИсаНоп$ Ва Бе допе. ТВе таш рагатаетгз оЁ агуше едитртеп угаз 
шуезбоаеа Бу сотрщег ехрегипеп 5. 

МешоЧ1са1| гесоттепдайоп Юг арр|уше оЁ еогейса! гезий аге дезспьеа. 

Вези! 5: ТВе Маш тези саПе4 оп айог$ оЁ Фе ага 15 ргезещед Бу зимезс 
Фтесйоп ш Фе Не!4 оЁ 4еуе!ортепЕ о {есВпо]оу оЁ агуте ади1лсКзап4 таепа|$ оп фе 
Базе оГаррИсаНоп «БоШи Тауег» уау. 

Сопеш$1оп: ТВе сотрайзоп шуезисайоп оЁ рагатеег$ оЁ агуше ргосез$ о п1015 
афилсКзап4 таепа! ш "БоШшс |ауег" оп тафетайса| по4е|$ ап4 1 геа! соп@ 1 юоп$ Угаз 
ргоу14е4. РтооЁ, а таетайса| то4е15 оГ деегтшеа фуре, Базе оп тафетайса| 
рВуз1с$ едпайоп$, у зайзКу |еуе] оф ргес1з1оп айгас@ сВагасег оЁ ргосез$ ап тау Бе 
тесоттепа Гог аррИсайоп. 
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